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Численное исследование сосуществования популяций в одной экологической нише* 


А. В. Будянский, М. Г. Кругликов, В. Г. Цибулин 


Описывается взаимодействие популяций хищников и жертв на пространственно неоднородном двумерном 
ареале. Модель записывается в виде системы нелинейных уравнений параболического типа для двух близко- 
родственных популяций хищников и двух популяций жертв, конкурирующих за общий ресурс. Показано, что 
при определённых соотношениях между параметрами и переменной по ареалу функции ресурса, модель при- 
надлежит к классу косимметричных динамических систем. В этом случае возникает непрерывное семейство 
стационарных распределений сосуществующих популяций. Вычислительный эксперимент основан на методе 
прямых и схеме смещённых сеток. Для аппроксимации по пространственным переменным задачи на прямо- 
угольном ареале используется метод баланса. Представлены результаты, демонстрирующие возможности мо- 
дели для описания формирования стационарных распределений популяций. Изучено формирование биологи- 
ческих структур при неоднородности параметров роста, проанализированы условия сосуществования близко- 
родственных видов. 

Ключевые слова: популяционная динамика, метод прямых, нелинейные параболические уравнения, косим- 
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Введение. Изменение и сокращение среды обитания биологических популяций в современном 
мире вызывает миграцию животных и приводит к смещению экологических равновесий. В процессе 
жизнедеятельности биологических видов образуются зоны совместного обитания (сосуществова- 
ния) популяций, причём размер этих зон зависит от интенсивности миграции, а процесс формиро- 
вания может быть достаточно медленным. Задача сохранения природного многообразия требует 
развития методов анализа динамики популяций хищников и их жертв [1]. 

Имеются различные точки зрения на возможность присутствия нескольких близкородствен- 
ных видов в одной экологической нише [1-4]. Принцип Гаузе [2] утверждает, что устойчивое сосу- 
ществование двух популяций невозможно, если рост ограничен одним жизненно важным ресурсом. 
В то же время, известны примеры [3], когда в одной экологической нише обитают несколько близ- 
кородственных популяций. В [4] показано, что при моделировании с учётом нелинейности мигра- 
ционных потоков возможна конкуренция биологических видов без вытеснения менее приспособ- 
ленной популяции (обобщение принципа Гаузе). 

В данной работе рассматривается модель, описывающая взаимодействие популяций хищ- 
ников и жертв. Для решения системы нелинейных уравнений в частных производных применяется 
метод конечных разностей. Целью работы является моделирование сосуществования видов и опре- 
деление параметров системы, при которых модель принадлежит классу косимметричных динами- 
ческих систем [5, 6] с непрерывным семейством стационарных распределений [7]. 

Модель динамики конкурирующих популяций. Рассматривается модель взаимодействия двух 
популяций жертв и двух видов хищников, представляющая собой систему параболических уравне- 
ний [8, 9]. Для описания динамики близкородственных популяций аналогично [10] используется 


единая функция обобщённого ресурса (ёмкости среды) р(х, И). Изменение плотности популяций 
жертв определяется логистическим законом (параметры роста |,, | ) и убылью из-за присутствия 


хищников (слагаемые с коэффициентами /,, /=1...4 ) [1]. 





* Работа выполнена при финансовой поддержке гранта РФФИ № 14-01-00470. 
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Прирост плотности популяций хищников за счёт потребления жертв даётся слагаемыми с 
коэффициентами |,, / =3...6 , аестественная смертность — слагаемыми с коэффициентами /., /, 
[11]. В (1) матрицы второго порядка к, состоят из диффузионных коэффициентов. 

Рассматривается ареал в виде прямоугольника © = |0,а ]х [0, В |, на границах которого ста- 
вились условия отсутствия потоков [12]: 

Уи, (0, у,Ё) =Уи, (а,у,Ё)=Уи, (х,0,Ё) =Уи, (х,Б,Ё)=0, Уи, (х,0,Ё) =Уи, (х,Б,Ё)=0, 
Уи, (0, у,Ё) =Уи/, (а,у,Ё) =0, Уи, (х,0,Ё) =Уи/,, (х,6,Ё) =0,1=1,2. 
Система дополняется начальными распределениями для плотностей популяций: 
0 0 г 
и, (х,у,0)=и; (х,у,0), и, (х,у,0)=и; (х,у,0), 1=12. (3) 

Рассматриваемая задача при определённых комбинациях параметров относится к классу 
косимметричных систем [5, 6], для которых возможно возникновение непрерывных семейств ста- 
ционарных распределений популяций [13]. 

Анализ показывает, что система (1) — (3) обладает косимметрией вида 

7. 
1 = (6.45, >И, › бЗИИЬ › ИИ, , (4) 
когда вещественные параметры с, связаны соотношениями, получающимися из условия ортого- 


(2) 


нальности вектора / правой части системы (1) — (2). Рассматриваемая система обладает косиммет- 
рией [14] при выполнении соотношений: 


К,» = Кб: и 46> = Моб,» 6» =—/ 61, 56» = 461, 


К.б4 =-К4бз, Изба = -М5бз, Мабд = Моб, 564 = —Кбз. 
Так как косимметрия [10] определяется с точностью до постоянного множителя, то далее 
СЕК —К 


эл” бз = Ида. 


(5) 


Численный метод. Для решения задачи (1) — (3) применяется метод прямых с дискретизацией 
на основе смещённых сеток. По переменным хи увводятся равномерные сетки: х, = 1, , У. = $й,, 
г=-1,0,...,п, +1, $=-1,0,...,, +1, №, =а/п,, №, =В/п, . Через и, „ (и, „) обозначается зна- 
чение плотности распределения популяции и, (и,) в узле [ХЕ У: Для вычисления потоков вво- 
дятся вспомогательные сетки: х,./› =-/, /2+Т,, г=1...., Пу, Ур =-Й, [2+51,, $=\...,П,. 


Компоненты Уи, (Уи/,) и Уи, (Уи) определяются в узлах (х я и [ху И соответ- 
г 5+ 


ственно. 
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При аппроксимации уравнений (1) по пространственным переменным вводятся разностные 
операторы первого порядка на двухточечных шаблонах: 


0215 - М, ; 051 — №, 5 
(9,6) „„, ы р ‚ (4,6), = Г . 


х У 
В результате получается следующая система обыкновенных дифференциальных уравнений: 
ее а 


РЕ. (6) 


С 
Дискретные аналоги краевых условий (2) дают соотношения для величин в законтурных 


узлах: 


Ил; Маз Ив 2И/,1;, $ = 0,..., п, 
(7) 
Иа ЕИ, аи И, АИ, Г=0,....П,, Г=12. 
Из (3) следуют начальные условия для (6), (7): 
Жи аи 
1=1,2, Г=0,...,П,, $ =0,...,П,. (8) 
Система (6) — (8) записывается в виде: 
ЕР О 
Ро (ИМ м, р 

и = ПАТИ "Или 22 1-1 пул, | ' (9) 


ии, сни ОНИ р } Г=1,2, 
,п,П, 


1,11 
И 


где У — вектор начальных данных. 


Для интегрирования по времени задачи Коши (9) применяется метод Рунге — Кутты чет- 
вёртого порядка. 

Косимметрией системы (9) является дискретный аналог векторного поля /[, получаемый из 
(4) — (5) в результате дискретизации. Поскольку нулевые решения задачи (9) (вектор У =0) ан- 


нулируют косимметрию, то любое ненулевое стационарное решение У, (то есть Е (У.) =0) не об- 
нуляет косимметрию и, таким образом, принадлежит однопараметрическому семейству равнове- 
сий. При этом в спектре устойчивости равновесия У. имеется нулевое значение, которое соответ- 


ствует нейтральному направлению вдоль семейства. Если остальные спектральные величины ле- 
жат в левой полуплоскости, то равновесие У, устойчиво. Это отвечает устойчивости в трансвер- 


сальном к семейству многообразии [5]. 
Численные результаты. Для прямоугольной области © = [0,2 | > [0,1] представлены результаты 


по исследованию стационарных распределений популяций жертв и хищников на прямоугольном 
ареале. Функция обобщённого ресурса р(х,у) соответствует наличию на ареале двух благопри- 


ятных зон (рис. 1). Расчёты проводились до выхода на устойчивые стационарные распределения. 
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В вычислениях использовались следующие па- 
раметры диффузии и роста-убыли: 


Ей 0 Е 0 
0 0,03 0 0,04 


Рис. 1. Функция обобщённого ресурса р (х,у) ры Ё 06 0 } в [г 08 0 
0 0,06 0 0,08 
ш =З, | =4, №: =З, 4. =1,2, Н; =4, № =1,6; 
ЕО Е 1, 2 
В этом случае система (1) — (3) обладает косимметрией (4) — (5) и возникает непрерывное семей- 
ство стационарных распределений. Для аппроксимации кривой семейства было вычислено более 
сотни точек, каждая из которых отвечает решению системы (9). Расчёт семейства проводился ме- 
тодом, описанным в [13]. В спектре полученных стационарных распределений имеются практиче- 
ски нулевые собственные значении (в =10`°). Это означает, что данные решения входят в нетри- 
виальное семейство стационарных распределений. При указанных параметрах семейство состоит 
из устойчивых решений, одно из них приведено на рис. 2. В силу однородности диффузии, профили 
распределений популяций жертв повторяют функцию обобщённого ресурса (см. рис. 1). При этом 
суммарная плотность жертв меньше значений функции р(х, у) в точках максимума и больше — в 











точках минимума. 

Если условие косимметрии нарушено, то для постоянных по пространству параметров роста 
происходит вытеснение наименее приспособленной популяции и на ареале остаётся только один 
вид. Например, если |. К, ., > Н.К, „ , ТО из двух популяций жертв выживает только и. 









0.2 < 
О 
а 





Рис. 2. Входящее в семейство стационарное распределение 
сосуществующих популяций жертв (а, Ь) и хищников (с, @) 


Сосуществование популяций без выполнения условий косимметрии возможно также при 
неравномерных по ареалу параметрах роста и убыли. Далее представлены результаты для пере- 


„ 21Х Н 
менного параметра роста популяции жертвы и, , задаваемого в виде р, = м», +н,, яп—— . В расчё- 
а 


тах было фиксировано значение п» =4. 
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1 2 3 4 5 6 м 
Рис. 3. Карта параметров, отвечающих сосуществованию (ПТ) и выживанию 
жертвы и, (Т) или и, (П), точка А — существование непрерывного семейства решений 


На рис. 3 дана карта режимов — плоскость параметров роста |, и модуляции |›,, которая 
состоит из трёх областей: две области соответствуют выживанию популяции и, или и, , а третья 


является областью сосуществования видов. Точка А(3,0), к которой «стягиваются» области, отве- 
чает случаю косимметрии и существованию семейства решений. 





Рис. 4. Стационарные распределения сосуществующих популяций в случае неравномерного по пространству параметра 
роста |, : их =4, в», =4 (столбец В); и, =4, м», =-3 (столбец С) 


На рис. 4 приведены распределения популяций для параметров |, =2,[,, =4 и 
ы =4, 5, =-2 (точки Ви Сна рис. 3). В этих случаях наблюдается устойчивое сосуществование 
популяций жертв и хищников на ареале, но в силу неоднородности параметра роста |, происходит 
размежевание популяций жертв. При |,, >0 (точка В) популяция и, доминирует в подобласти 
х<а/2, а популяция и, — в остальной части ареала. Плотности распределений популяций хищ- 
ников в общем напоминают профиль популяции и, ‚ что связано с большими коэффициентами роста 
за счёт популяции и, . При р,, <0 (точка С) картина изменяется: теперь и, доминирует в подоб- 


ласти х>а/2. 
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В случае неравномерности параметров роста хищников по ареалу, также реализуются сце- 
нарии сосуществования и вытеснения видов. Для неоднородного параметра роста 


. 21х 
М5 = Му + В, Яп-— при фиксированном значении п.’ = 4 на рис. 5 приведена карта режимов на 
а 


плоскости параметра |, и коэффициента модуляции р... Карта состоит из трёх областей: области 


совместного существования хищников и областей выживания одного из хищников. Границы области 
стягиваются к точке, удовлетворяющей условию косимметрии (р., =0, м. =3). 


Заключение. Найдены условия на параметры системы, при которых модель принадлежит классу 
косимметричных динамических систем [10, 11] и имеется непрерывное семейство стационарных 
распределений [14]. Предложена модель для описания взаимодействия популяций хищников и 
жертв на пространственно неоднородном двумерном ареале. Модель представлена в виде системы 
нелинейных уравнений параболического типа для двух близкородственных популяций хищников и 
двух популяций жертв, конкурирующих за общий ресурс. 
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Рис. 5. Карта параметров, отвечающих сосуществованию (11) и 
выживанию одного их хищников и/, (Ги и, (П) 


Из 


Показано, что при определённых соотношениях между параметрами и переменной по аре- 
алу функции ресурса, задача принадлежит классу косимметричных динамических систем. В этом 
случае возникает непрерывное семейство стационарных распределений сосуществующих популя- 
ций. Вычислительный эксперимент основан на методе прямых и схеме смещённых сеток, аппрок- 
симации построены с использованием метода баланса. Результаты демонстрируют возможности 
модели для описания формирования стационарных распределений популяций. Проанализированы 
условия сосуществования близкородственных видов, изучено формирование биологических струк- 
тур при неоднородности параметров роста. 

Библиографический список 

1. Миггау, 3. О. Мабетабса! Вюоду Ц. браца! тоде!5 апа Воюте са! АррйсаНопз / 3. В. Миг- 
гау. — 5рипдег—\Мейад, 2003. — 1082 р. 

2. Гаузе, Г. Ф. Борьба за существование / Г. Ф. Гаузе. — Ижевск : Ин-т компьютерных ис- 
следований, 2002. — 234 с. 

3. Бигон, М. Экология. Особи, популяции и сообщества / М. Бигон, Дж. Харпер, 
К. Таунсенд. — Москва : Мир, 1989. — 1144 с. 

4. Белотелов, Н. В. Популяционные модели с нелинейной диффузией / Н. В. Белотелов, 
А. И. Лобанов // Математическое моделирование. — 1997. — Т. 9, № 12. — С. 43-56. 

5. Юдович, В. И. Косимметрия, вырождение решений операторных уравнений, возникно- 
вение фильтрационной конвекции / В. И. Юдович // Математические заметки. — 1991. — Т. 49, 
№ 5. — С. 142-148. 

6. Учао\си, \. Т. бесопдагу суще оЁ едина т а эзует мА созутте{у, К5 сгеабоп Бу 
ЫгРгсайоп апа итроз$ЬИКу ог зуттей@с {геайтепЕ о! К. СПаоз, 1995, \о1. 5, по. 2, рр. 402—411. 


33 


Математические и естественные науки 








7. бомогикРт, \/. Сотриег ехрейтепт 5 мИИ созуттей1с тоаек. 7. Апдеми. Мат. Месн, 1996, 
уо|. 76, рр. 559-562. 

8. Вапеде, М., Регоу$К, $5. 5еР-огдатзе4 зраЧа!| раЦегп$ апа сНао$ т а гаНо-дерепаеч ргед- 
аюг-ргеу зу$ет. 2. ТВеог. Вю|., 2011, \а.. 4, рр. 37-53. 

9. Хие, Е. РаКегп ГогтаНоп т а ргедатог-ргеу тоде! млн раба! еНесе. РНууса А: ЗбаНИса! 
теспап!с$ апа #5 аррИсаНоп$, 2012, у. 391, рр. 5987—5996. 

10. Будянский, А. В. Моделирование пространственно-временной миграции близкородствен- 
ных популяций / А. В. Будянский, В. Г. Цибулин // Компьютерные исследования и моделирование. — 
2011. — Т. 3, № 4. — С. 477—88. 

11. Мишугова, Г. В. Моделирование процесса загрязнения атмосферы / Г. В. Мишугова // 
Вестник Дон. гос. техн. ун-та. — 2012. — № 8 (69). — С. 12-17. 

12. Заковортный, В. Л. Моделирование эволюции динамической системы, взаимодействую- 
щей со средой / В. Л. Заковоротный, Фам Дин Тунг // Вестник Дон. гос. техн. ун-та. — 2006. — Т. 6, 
№ 3 (30). — С. 184-200. 

13. КоуаЕетма, Е. $., РИзсНтиЕй, К., ТзуБийт, \. С. Вупаптс$ ог попИтеаг рагабойс едцаНоп$ МН 
созуттейу. Сотрщег Ацега т 5аепИЯс Сотрийпа, САЗС, 2007, рр. 265—274. 

14. РАзситиКИ, К., Комаема, Е. $., ТзуБийт, У. С. Ратйу оЁ едина т а рорМ4аНоп Кте{с$ 
поЧе!| апа №5 со!арзе. МопИпеаг Апа!у$5: Кеа! \М/ойа АррйсаНопз, 2011, уа.. 12, рр. 145-155. 


Материал поступил в редакцию 01.11.2013. 


ВеГегепсе$ 

1. Миггау, 2. О. МалетаЧса! Воду П. 5раНа! поде|5 апа ВотеФ!са! АррИсаНопз. 5рипдег— 
\епад, 2003, 1082 р. 

2. Саицге, С. Е. ВогЬа га зизНспе$Уоуатуе. [5гидое Гог ех$епсе.] 12НемзК: Тпьи 
КотрущегпукВ 53 едомапу, 2002, 234 р. (м Визчап). 

3. Ведоп, М., Нагрег, 2., Томтзепа, С. Екоодуа. ОзоЫ, рору!уа{ 1 зооб$ПсНезма. [Есооду. 
пам@чца, рорЧаНоп$ апа соттипфез.] Мозсоми: Ми, 1989, 1144 р. (т Визчап). 

4. Васею\м, М. \., ГоБапох, А. Т. Рориуаоппууе тодей $ пе!теупоу ЧЁитуеу. [РормаНоп 
тоде!5 мИН поп-!пеаг ЧЁизюп.] Маета{спезКоуе тодегоуапуе, 1997, \уо|. 9, по. 12, рр. 43—56 (м 
Киап). 

5. Учаомси, \. 1. Козттеуа, уугогПаетуе гезпету орегаогпуКН игаупету, уо7тппоуетуе 
Яга юоппоу КопуеКй. [Созутте!у, дедепегайоп оЁ орегаюг едиаНоп зоиЧоп$, опзеё оЁ ЯЁгаНоп 
сопуесЧоп.] Маетайсйе$Куе хате 4, 1991, уо/. 49, по. 5, рр. 142-148 (т Визчап). 

6. Уиаомсй, \. Т. бесопдагу суще оЁ едийБпа т а зу$ет мВ созуттейу, К$ сгеабоп Бу 
ЫгРигсайоп апа итроз$ИКу ог зуттей@с {геайтепЕ о! К. СПаоз, 1995, \о1. 5, по. 2, рр. 402—411. 

7. СомогикИт, \. Сотрщег ехрейтепт5 уИ созутте с тоде{5. 2. Апдеми. Ма{1. Месп, 1996, 
уо|. 76, рр. 559-562. 

8. Вапеде, М., Регоу$К, 5. 5е-огдат5е4 зраНа! раЦегпз апа сПао$ п а гаЧо-дерепаеч ргед- 
аюг-ргеу эу%ет. 2. ТВеог. Вю|., 2011, ма. 4, рр. 37-53. 

9. Хие, Е. РаЦегт ГогтаНоп шт а ргедаюг-ргеу по! мН зраЧа! еКесе. Рпузса А: ЗаНзНса| 
теспап!с$ апа #5 аррИсаНоп$, 2012, у. 391, рр. 5987—5996. 

10.ВидуапзКу, А. \М., ТзуБи!т, \. С. Модейоуатуе ргобгап$Меппо-угетеппоу плюгай Ып- 
Когод$\/еппукй рорцуау. [Модей!пд оЁ зраНа!-{етрога! птогабоп оЁ ‹озеу гаае рор\аНопс.] 
Котрущегпууе 15$едоуатуа 1 поде!йоуатуе, 2011, ума. 3, по. 4, рр. 477—488 (т Визчап). 

11.Мпидома, С. \. Аг соаттайоп ргосез$ эитиа{оп. Мезп оЁ О$ТЦ, 2012, по. 8 (69), 
рр. 12-17 (м ВизЗап). 


34 


Вестник ДГТУ. 2014. Т. 14, №2 (77) 








12.Гакомогопу, \. Е., РВат ОВ Типд. Модейгоуатуе емоуцй АтаписпезКоу у$ету, ута!- 
тодеу$иучзЙспеу 50 эгедоу. [5итМаЧоп оЁ еуо!иНоп оЁ дупаптса! зу$ет п\{егас/та мий теит.] 
\Мезп оЁ О$ТЦ, 2006, \о1. 6, по. 3 (30), рр. 184-200 (т Виззап). 

13.Комае\а, Е. $., РИ5ситтиКИ, К., ТзуБийт, \. С. Вупапс$ оЁ поп!пеаг рагабойс едцаНоп$ мИП 
созуттейу. Сотрщег Ацебга т 5аепИЯс Сотрийпа, САЗС, 2007, рр. 265—274. 

14.Рй5сбтиёН, К., Коуаема, Е. $., ТзуБийт, У. С. Ратйу ог едийБйа т а рормаНоп Кте{с$ 
поЧе!| апа №5 со!арзе. МопИпеаг Апа!у$5: Кеа! \/ойа АррйсаНопз, 2011, уа.. 12, рр. 145-155. 


МОМЕКТСАЕ 5ТООУ ОЕ СОЕХТЗТЕМСЕ ОЕ РОРИЕАТТОМ$ 1М АМ ЕМУТКОММЕМТАЕ МТСНЕ* 
А. \. Видуап$Ку, М. С. КгидИКох, У. 6. Тзубийт 


Тре ргедайог-ргеу Гпёегасйоп$ оп {Ве рава! ВЕегодепеоиз #вио-айтепзюпа! агеа аге аезспреа. Тпе тодЕ!//5 ииеп 
аз а 5у$ет оЁ поп/теаг рагабо!с едиайоп$ Гог о сое геаёеа ргедаюг рорш/айоп$ апа #/ю ргеу роршайоп5 
сотрейпд Юг {те депега! гезоигсе. [& 15 зпоит {паЁ ипаег сейат ге/авоп$Шрз Бейееп те рагатеег5 апа {пе 
уама Ме паёига/! ВабаЕ гезоигсе ипсНоп5, те тодЕ! Беопдс №0 {пе са5$ ое созуттей\т аупатиа! зу$етэ. п 
{1/5 сазе, {пеге [5 а сопйпиоиз гатйу оЁ $авопагу @5БинНоп$ о пе соежептЕ рори/айоп5. ТВе зитшайоп ехрег- 
/тепЕ 65 Базе оп {Пе теёоа о! га/дйЕ те, апа оп {те °срете о! Яаоддегеа ата. Тре Ба/апсе те!тоа 5 изед 
Юг Пе арргохмтайоп т зрайа! иапаМЕ$ ОЕ {Пе Еа5К оп а гецапдий! агеа. ТПе гезш!5 звои/пд {пе тоде! сара И#е5 
Юг аезстртд {те юЮгтайоп о! Не роршайоп чаНопагу @5трБиНопз аге ргезептеа. Тйе юптайоп о! пе Бо/оба! 
$гисигев 15 зи еа ипаЕег {пе дгоиЛй рагатеег пЕегодепейу; {те сопа!юопз Юг те соемепсе оЁ сое гевеа 
[урез аге апа\угеа, 

Кеуигога5: рори/айоп дупат/$, теЁпоа оЁ а/о ЙЕ [тез, поп/пеаг рагаро/!с едиаНоп5, созуттебу. 





* Тпе гезеагси 15 допе ми {Пе Ипапса! зирро[ гот ВЕЕТ (дгапе по. 14-01-00470). 
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